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Abstract: We consider the high-frequency Helmholtz equation with a given 
source term, and a small absorption parameter a > 0. The high-frequency (or: 
semi-classical) parameter is £ > 0. We let £ and a go to zero simultaneously. We 
assume that the zero energy is non-trapping for the underlying classical flow. 
We also assume that the classical trajectories starting from the origin satisfy a 
transversality condition, a generic assumption. 
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Under these assumptions, we prove that the solution radiates in the outgo¬ 
ing direction, uniformly in e. In particular, the function u®, when conveniently 
rescaled at the scale e close to the origin, is shown to converge towards the out¬ 
going solution of the Helmholtz equation, with coefficients frozen at the origin. 
This provides a uniform version (in e) of the limiting absorption principle. 

Writing the resolvent of the Helmholtz equation as the integral in time of 
the associated semi-classical Schrodinger propagator, our analysis relies on the 
following tools: (i) For very large times, we prove and use a uniform version 
of the Egorov Theorem to estimate the time integral; (ii) for moderate times, 
we prove a uniform dispersive estimate that relies on a wave-packet approach, 
together with the above mentioned transversality condition; (iii) for small times, 
we prove that the semi-classical Schrodinger operator with variable coefficients 
has the same dispersive properties as in the constant coefficients case, uniformly 
in e. 

2000 Mathematics Subject Classification number: Primary 35Q40; Sec¬ 
ondary 35J10, 81Q20. 


1 Introduction 

In this article, we study the asymptotics e ^ 0'*' in the following scaled Helm¬ 
holtz equation, with unknown rc®, 

ie Oe w^(x) -I- ^Aa;'w^{x) + n?{ex)w‘^ {x) = S{x). (1.1) 

In this scaling, the absorption parameter Og > 0 is small, i.e. 

Oe ^ O'*' as £ ^ 0. 

The limiting case = O’*" is actually allowed in our analysis. Also, the index 
of refraction n?(ex) is almost constant, 

n^{ex) « n^(0). 

The competition between these two effects is the key difficulty of the present 
work. 

In all our analysis, the variable x belongs to for some d > 3. The index 
of refraction (a:) is assumed to be given, smooth and non-negative^ 

Vx G n^{x) > 0, and n^{x) G (1.2) 

It is also supposed that n^(x) goes to a constant at infinity, 

n^{x) = + O (^{x)~^) as x ^ oo, (1.3) 

for some, possibly small, exponant p > 0^. In the language of Schrodinger 
operators, this means that the potential — n'^{x) is assumed to be either 

^Our analysis is easily extended to the case where the refraction index is a function that 
changes sign. The only really important assumption on the sign of n is > 0, see Proposition 
Otherwise, all the arguments given in this paper are easily adapted when changes 

sign, the analysis being actually simpler when ri^{x) has the wrong sign because contribution 
of terms involving xsi^e) vanishes in that case (see below for the notations). 

^Here and below we use the standard notation {x) := (1 + . 
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short-range or long range. Finally, the source term in uses a function S{x) 
that is taken sufficiently smooth and decays fast enough at infinity. We refer 
to the sequel for the very assumptions we need on the refraction index ri^{x), 
together with the source S (see the statement of the main Theorem below). 

Upon the L^-unitary rescaling 

w^{x) = e’^/‘^u^{ex), 

the study of cm is naturally linked to the analysis of the high-frequency Helm¬ 
holtz equation, 

ieasU^x) + y^xU^x) + n'^{x)u^x) = , (1.4) 

where the source term S{x/e) now plays the role of a concentration profile at the 
scale e. In this picture, the difficulty now comes from the interaction between the 
oscillations induced by the source S{x/e), and the ones due to the semiclassical 
operator e^A/2 -|- ii?{x). We give below more complete motivations for looking 
at the asymptotics in dJ or 

The goal of this article is to prove that the solution to dJ converges 
(in the distributional sense) to the outgoing solution of the natural constant 
coefficient Helmholtz equation, i.e. 

lim , where is defined as the solution to 

e^O 

iO+r(;™‘(a;) -f iA,,r(;™*(a;) -b {0)w°'^Hx) = S (x). (1.5) 

In other words, 

^out _ g 

5^0+ V 2 J 

= i j exp ^^Aj,-b n^(0)^^ S'dt. (1.6) 

It is well-known that u>°“* can also be defined as the unique solution to (Aa;/2 -b 
n^(0))ui°'^‘ = S that satisfies the Sommerfeld radiation condition at infinity 

-^^■W^w°'^^{x)+in{0)w°'^^(x)=o(j^'^ as |a;| ^ oo. (1.7) 

The main geometric assumptions we need on the refraction index to ensure 
the validity of dll are twofolds. First, we need that the trajectories of the 
Hamiltonian ^^/2 — n^{x) at the zero energy are not trapped. This is a stan¬ 
dard assumption in this context. It somehow prevents accumulation of energy in 
bounded regions of space. Second, it turns out that the trajectories that really 
matter in our analysis, are those that start from the origin x = 0, with zero en¬ 
ergy ^^/2 = n^(0). In this perspective, we need that these trajectories satisfy a 
transversality condition: in essence, each such ray can self-intersect, but we 
require that the self-intersection is then “tranverse” (see assumption (H) page 
E3 i-e. dSli, , in section 0 below). This second assumption prevents 

accumulation of energy at the origin. 

We wish to emphasize that the statement dl is not obvious. In particular, 
if the transversality assumption (H) page|23is not fullfilled, our analysis shows 
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that dni) becomes false in general. We also refer to the end of this paper for 
“counterexamples”. 

The central difficulty is the following. On the one hand, the vanishing ab¬ 
sorption parameter in im leads to thinking that should satisfy the Som- 
merfeld radiation condition at infinity with the variable refraction index 
n'^{sx) (see lll.7ll l. Knowing that this roughly means 

that w® should behave like exp(i2“^/^noo|a^|)/|a;| at infinity in x (in dimension 
d = 3, say). On the other hand, the almost constant refraction index n^(£x) in 
inj leads to observe that naturally goes to a solution of the Helmholtz equa¬ 
tion with constant refraction index n?{0). Hoping that we may follow the 
absorption coefficient continuously along the limit e ^ 0 in n'^{ex), the state¬ 
ment C31) becomes natural, and w'^ should behave like exp(*2 ^/^n(0)|a;|)/|x| 
asymptotically. But, since n(0) ^ Uoa in general, the last two statements are 
contradictory ... As we see, the strong non-local effects induced by the Helmholtz 
equation make the key difficulty in following the continuous dependence of w'^ 
upon both the absorption parameter —*■ 0+ and on the index n'^{ex) —> n^(0). 

Let us now give some more detailed account on our motivations for looking 
at the asymptotics £ ^ 0 in CU. 

In [BCKP) . the high-frequency analysis of the Helmholtz equation with 
source term is performed. More precisely, the asymptotic behaviour as e —> 0 
of the following equation is studied^ 

ieasU^x) + ^/^^u’^{x) + n'^{x)u^x) = , (1.8) 

where the variable x belongs to for some d > 3, and the index of re¬ 
fraction n^(x) together with the concentration profile S(x) are as before (see 
|BCKP| 1. Later, the analysis of |B(1KP| was extended in |(1PH| to more gen¬ 
eral oscillating/concentrating source terms. The paper ITTFm studies indeed 
the high-frequency analysis e ^ 0 in 

£2 

ieaeU^{x) -I- —AxU^{x) +n'^{x)u^{x) = 

^( 2 /) da{y). (1.9) 

(See also |C^ for extensions - see for the case where n? has discontinu¬ 
ities). In H1.9|l . the function S again plays the role of a concentration profile like 
in but the concentration occurs this time around a smooth submanifold 

r C of dimension p instead of a point. On the more, the source term here 
includes additional oscillations through the (smooth) amplitude A and phase 
(j). In these notations da denotes the induced euclidean surface measure on 
the manifold P, and the rescaling exponant q depends on the dimension of P 
together with geometric considerations, see iHFm . 

Both Helmholtz equations CSJ and modelize the propagation of a 

high-frequency source wave in a medium with scaled, variable, refraction index 
in?{x)/s‘^. The scaling of the index imposes that the waves propagating in the 
medium naturally have wavelength e. On the other hand, the source in 

^note that we use here a slightly different scaling than the one used in |RCKP| . This a 
harmless modification that is due to mere convenience. 
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as well as is concentrating at the scale £, close to the origin, or close to 
the surface F. It thus carries oscillations at the typical wavelength e. One 
may think of an antenna concentrated close to a point or to a surface, and 
emmitting waves in the whole space. The important phenomenon that these 
linear equations include precisely lies in the resonant interaction between 
the high-frequency oscillations of the source, and the propagative modes of the 
medium dictated by the index v? je^. This makes one of the key difficulties of 
the analysis performed in |BOKP| and IHFrI . 

A Wigner approach is used in [BCIKP] and |( fPR] to treat the high-frequency 
asymptotics £ ^ 0. Up to a harmless rescaling, these papers establish that the 
Wigner transform of u^{x) satisfies, in the limit £ ^ 0, the stationnary 

transport equation 

0+/(a:, 0 + e ■ V,/(x, 0 + V,n2(x) • V^fix, 0 = Qix, 0, (1-10) 

where /(x,^) = lim/^(x,^) measures the energy carried by rays located at 
the point x in space, with frequency ^ G The limiting source term Q 

in firnii describes quantitatively the resonant interactions mentioned above. 
In the easier case of (11.811 . one has Q{x,^) = (5 (^^/2 — n^(0)) S{x) |S'(C)Pj 
meaning that the asymptotic source of energy is concentrated at the origin in 
X (this is the factor <5(x)), and it only carries resonant frequencies ^ above this 
point (due to 6 (^^/2 — n^(0))). A similar but more complicated value of Q is 
obtained in the case of In any circumstance, equation n~Tnii tells us that 
the energy brought by the source Q is propagated in the whole space through 
the transport operator ^ -I- Va;n^(x) • Vj naturally associated with the semi- 

classical operator —e'^Ax/2 — n'^{x). The term O'*"/ in n~Tnii specifies a radiation 
condition at infinity for /, that is the trace, as £ ^ 0 of the absorption coefficient 
«£ > 0 in Oil and Oil . It gives / as the outgoing solution 


r+oo 

f{x,0= Q(Af(s,x,^),S(s,x,^)) ds. 

Jo 


Here (A(s, x, ^), S(s, x, ^)) is the value at time s of the characteristic curve of 
^ • Va; -I- VxU^(x) • starting at point (x,^) of phase-space (see (II.UH below). 
Obtaining the radiation condition for / as the limiting effect of the absorp¬ 
tion coefficient in is actually the second main difficulty of the analysis 
performed in |BCKP| and EES. 

It turns out that the analysis performed in IBOKPI relies at some point 
on the asymptotic behaviour of the scaled wave function w®(x) = (ex') 

that measures the oscillation/concentration behaviour of close to the origin. 
Similarly, in ITWl one needs to rescale u® around any point y G F, setting 
Wy(x) := £‘^/^u®(y -I- ex) for any such y. We naturally have 

ieaeW^(x) + ^Axw‘^(x) + n^(£x)w®(x) = S (x), 

in the case of and a similar observation holds true in the case of 

Hence the natural rescaling leads to the analysis of the prototype equation 
ini. Under appropriate assumptions on n^(x) and S(x), it may be proved 
that w'^ , solution to lEU, is bounded in the weighted space L^((x)^+'^ dx), 
for any d > 0, uniformly in e. For a fixed value of £, such weighted estimates are 
consequences of the work by Agmon, Hormander, EHi The fact that these 
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bounds are uniform in e is a consequence of the recent (and optimal) estimates 
established by B. Perthame and L. Vega in IFTTI . EYl (where the weighted 
L? space are replaced by a more precise homogeneous Besov-like space). The 
results in lES] and IpH actually need a virial condition of the type 2'n?'(x)+x- 
\/xn^(x) > c > 0, an inequality that implies both our transversality assumption 
(H) page ESI and the non-trapping condition, i.e. the two hypothesis made 
in the present paper. We also refer to the work by N. Burq |Bu| . Gerard and 
Martinez |GM |. T. Jecko m , as well as Wang and Zhang |WZ| , for (not optimal) 
bounds in a similar spirit. Under the weaker assumptions we make in the present 
paper, a weaker bound may also be obtained as a consequence of our analysis. 
In any case, once w'^ is seen to be bounded, it naturally possesses a weak limit 
w — limrc'^ in the appropriate space. The limit w clearly satisfies in a weak 
sense the equation 

+ w{x) = S{x). (1-11) 

Unfortunately, equation does not specify w = limru^ in a unique way, 

and it has to be supplemented with a radiation condition at infinity. In view 
of the equation cu satisfied by w^, it has been conjectured in |BGKP| and 
(TWI that limw^ actually satisfies 

lim w'^ = w™*, 

where is the outgoing solution defined before. The present paper answers 
the conjecture formulated in these works. It also gives geometric conditions for 
the convergence lim w® = w°'^^ to hold. 

As a final remark, let us mention that our anaylsis is purely time-dependent. 
We wish to indicate that similar results than those in the present paper were re¬ 
cently and independently obtained by Wang and Zhang |WZ| using a stationary 
approach. Note that their analysis requires the stronger virial condition. 


Our main theorem is the following 

Main Theorem 

Let w’^ satisfy ieasw'^{x) + ^AxW^{x) +nf(ex)w^{x) = S{x), for some se¬ 
quence «£ > 0 such that Og —> O'*" as e —> 0. Assume that the source term S 
belongs to the Schwartz class ^(K^^). Suppose also that the index of refraction 
satisfies the following set of assumptions 

• (smoothness, decay). There exists an exponent p > 0, and a positive 
constant > 0 such that for any multi-index a G N^, there exists a 
constant Co, > 0 with 


9“ - n^) < Co (x) ^ '“L 


( 1 . 12 ) 


• (non-trapping condition). The trajectories associated with the Hamilto¬ 
nian ^ (2 — n^(x) are not trapped at the zero energy. In other words, any 
trajectory (V(t, cc, ^), S(t, a;, ^)) solution to 


d 

—X{t,x,f) = E{t,x,f,), 
^E{t,x,0 = {X{t,x,f)), 
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V(0,a;,^) = X, 
S(0,x,C) = 


(1.13) 
















with initial datum (x,^) such that — n^(x) = 0, is assumed to satisfy 


|X(i, a;,^)| ^ oo, as |i| ^ oo. 

• (tranversality condition). The tranvsersality condition (H) nn.ne \SfA fsee 
also and J 7. 2)^ ) on the trajectories starting from the origin x = 0, 

with zero energy /2 = n^(0), is satisfied. 

Then, we do have the following convergence, weakly, when tested against any 
function 4> € S(R‘^), 


First remark 

Still referring to (H) page|2Sl or 17.21111 . 17.2411 1 for the precise statements, we 
readily indicate that the transversality assumption (H) essentially requires that 
the set 

{{r,, t) G R2dx]0, oo[ s.t. X{t, 0, C) = 0, S(t, 0, ^ = V, eV2 = ^'(0)} 

is a smooth submanifold of having a codimension > d + 2, a generic 

asssumption. In other words, zero energy trajectories issued from the origin 
and passing several times through the origin x = 0 should be “rare”. □ 

Second remark 

As we already mentionned, it is easily proved that the virial condition 2n^(x) + 
X ■ Va:n^(x) > c > 0 implies both the non-trapping and the transversality con¬ 
ditions. This observation relies on the identities dt {X{t,x,f,)‘^/2) = X{t,x,f,) ■ 
E{t,x,^) and dt{X{t,x,J,) ■E{t,x,^)) = [2n^{x) + x-X a:n^{x)\\,,=x(t,x,0 > 
c > 0, where (X{t,x,ff),'E.{t,x,Sf) is any trajectory with zero energy (see section 
El for computations in this spirit). 

In fact, the virial condition implies even more, namely that trajectories issued 
from the origin with zero energy never come back to the origin. In other words, 
the set involved in assumption (H) page 03 is simply void, and (H) is trivially 
true under the virial condition. As the reader may easily check, such a situation 
allows to considerably simplify the proof we give here: the tools developped in 
sections EElISlil are actually enough to make the complete analysis, and one 
does not need to go into the detailed computations of section [3 in that case. 

Last, the above Theorem asserts the convergence of w^: note in passing 
that even the weak boundedness of w® under the sole above assumptions (i.e. 
without the virial condition) is not a known result. □ 

The above theorem is not only a local convergence result, valid for test 
functions (p G S. Indeed, by density of smooth functions in weighted spaces, 
it readily implies the following immediate corollary. It states that, provided rc® 
is bounded in the natural weighted space, the convergence also holds weakly 
in this space. In other words, the convergence also holds globally. 

Immediate corollary 

With the notations of the main Theorem, assume that the source term S above 
satisfies the weaker decay property 

\\S\\B-.= J22^/^S\\r,^c,)<oo, ( 1 . 14 ) 

i6Z 
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where Cj denotes the annulus {2^ < |a;| < 2-^+^} in Suppose also that the 
index of refraction satisfies the smoothness condition of the main Theorem, with 
the non-trapping and transversality assumptions replaced by the stronger 


• (virial-like condition) 


2E 




sup 

xGCj 


[x ■ Vn^(a;)) 
n‘^{x) 


< 1 . 


(1.15) 


Then, we do have the convergence weakly, when tested against any 

function (j) such that ||(^||b < oo. 

Under the simpler virial condition 2n^{x) + x ■ n^{x) > c > 0, a similar result 
holds with the space B replaced by the more usual weighted space (^{x^^^dx) 
((5 > 0 arbitrary). Here, we give a version where the decay 11.1411 assumed on 
the source S is the optimal one, and the above weak convergence holds in the 
optimal space. 

It is well known that the resolvent of the Helmholtz operator maps the 
weighted space to (^{x)~^~^dx^ for any d > 0 ( |Ag| , p], 

|(lMp . Agmon and Hormander na gave an optimal version in the constant 
coefficients case: the resolvent of the Helmholtz operator sends the weighted Lf 
space B defined in (tmii to the dual weighted space B* defined by 


Is* := sup2 


(1.16) 


For non-constant coefficients, that are non-compact perturbations of constants, 
Perthame and Vega in imi and |PV2| established the optimal estimate in B- 
B* under assumption n~T^ . In our perspective, the assumption 11.1511 is of 
technical nature, and it may be replaced by any assumption ensuring that the 
solution to cu satisfies the uniform bound 

Ww’^Wb- < C'd,„2 \\S\\b, (1-17) 

for some universal constant „2 that only depends on the dimension d > 3 
and the index n^. 

Proof of the immediate Corollary 

Under the virial-like assumption 11.1511 . it has been established in Em that 
estimate (im holds true. Hence, by density of the Schwartz class in the space 
B, one readily reduces the problem to the case when the source S and the test 
function (j> belong to 5(R'^). The Main Theorem now allows to conclude. □ 

Needless to say, the central assumptions needed for the theorem are the non¬ 
trapping condition together with the transversality condition. Comments are 
given below on the very meaning of the transversality condition (H) page 11151 
(i.e. rr^ . (TT^ l. to which we refer. 

To state the result very briefly, the heart of our proof lies in proving that un¬ 
der the above assumptions, the propagator exp {ie~^t {—e^/S.x/2 — n^{x))), or 
its rescaled value exp (it (—Aa;/2 — nf{ex))), satisfy “similar” dispersive prop¬ 
erties as the free Schrodinger operator exp {it (—A2,/2 — n^(0))), uniformly in e. 
This in turn is proved upon distinguishing between small times, moderate times, 
and very large times, each case leading to the use of different arguments and 
techniques. 









The remainder part of this paper is devoted to the proof of the main The¬ 
orem. The proof being long and using many different tools, we first draw in 
section |21 an outline of the proof, giving the main ideas and tools. We also 
define the relevant mathematical objects to be used throughout the paper. The 
proof itself is performed in the next sections EltolHl Examples and counterex¬ 
amples to the Theorem are also proposed in the last section |21 

The main intermediate results are proposition^ proposition|21 proposition^ 
together with the more difficult proposition 0] (that needs an Egorov Theorem 
for large times stated in Lemma The key (and most difficult) result is 
proposition [3 The latter uses the tranversality condition mentioned before. 


2 Preliminary Analysis: outline of the proof of 
the Main Theorem 

2.1 Outline of the proof 

Let w® be the solution to isa^w^ + ^ Aw® -|- n^(ea;)w® = S (x), with S G 5(]R‘^). 
According to the statement of our main Theorem, we wish to study the asymp¬ 
totic behaviour of w® as e —> 0, in a weak sense. Taking a test function 
4>{x) G 5(IR‘^), and defining the duality product 

(w®,0):= / {x)(j){x) dx, 

JR'^ 

we want to prove the convergence 

(w®, (^) ^ (w™*, (()) as e ^ 0. 

where the outgoing solution of the (constant coefficient) Helmholtz equation 
is defined in before. 

First step: preliminary reduction - the time dependent approach 

In order to prove the weak convergence we define the 

rescaled function 

(f) ■ (2-1) 

It satisfies ieae'«®-|-e^/2 AM®-|-n^(a;)u® = 1/e'^^‘^S {x/e) =: Se{x), where for any 
function f{x) we use the short-hand notation 

f.M = (j) ^ 

Using now the function u® instead of w®, we observe the equality 

(w®,0) = (M®,()ie). (2.2) 

This transforms the original problem into the question of computing the semi- 
classical limit e ^ 0 in the equation satished by u®. One sees in that this 
limit needs to be computed at the semiclassical scale (i.e. when tested upon a 
smooth, concentrated function ipg). 

In order to do so, we compute m® in terms of the semiclassical resolvent 
(teoE -I- (e^/2)A-|-n^(x)) . It is the integral over the whole time interval 
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[0, +oo[ of the propagator of the Schrodinger operator associated with £^A/2 - 
In other words we write 


= [ lea^ +—A + n'^{x) 


S, 


= i 


J exp +—A + n^{x)^'^ Sg dt. (2-3) 


Now, defining the semi-classical propagator 


t / e 


Ue{t) := exp A -h n (x) j j = exp , 

associated with the semi-classical Schrodinger operator 


(2.4) 


H,:=--A-n^{xl 


we arrive at the final formula 




^+oo 


{US)SeAe)dt. 


(2.5) 


( 2 . 6 ) 


Our strategy is to pass to the limit in this very integral. 


Second step: passing to the limit in the time integral (|2:6|l 

In order to pass to the limit e ^ 0 in (EH), we need to analyze the contribu¬ 
tions of various time scales in the corresponding time integral. More precisely, 
we choose for the whole subsequent analysis two (large) cutoff parameters in 
time, denoted by Tq and Ti, and we analyze the contributions to the time 
integral (EH that are due to the three regions 

0 < t < To e, Tq e < t < Ti, and t > Ti. 

We also choose a (small) exponent k > 0, and we occasionally treat separately 
the contributions of very large times 

t > e””. 


Associated with these truncations, we take once and for all a smooth cutoff 
function y defined on R, such that 

Xiz) = l when \z\<l/2, x(z) = 0 when |z| > 1, 

x{z) > 0 for any z. (2.7) 

To be complete, there remains to finally choose a (small) cutoff parameter in 
energy 5 > 0. Accordingly we distinguish in the scalar product {Ue{t)Se, (jie) 
between energies close to (or far from) the zero energy, which is critical for our 
problem. In other words, we set the self-adjoint operator 

(iJe) ■=x(^^ ■ 
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This object is perfectly well defined using standard functional calculus for self- 
adjoint operators. We decompose 

mt)S„cj,,) = (uS)xs{He)S,,<i,,) + (uS) (1 - Xi) 

Following the above described decomposition of times and energies, we study 
each of the subsequent terms: 

• The contribution of small times is 

1 e f t \ 

£ Jo ^ V^J 

We prove in section 13 that this term actually gives the dominant contribution in 
(EH), provided the cutoff parameter Tq is taken large enough. This (easy) anal¬ 
ysis essentially boils down to manipulations on the time dependent Schrodinger 
operator idt + Aa;/2 -|- ri^{ex), for finite times t of the order t ^ Tq at most. 

• The contribution of moderate and large times, away from the zero 
energy, is 

J (1 - x) “ ^' 5 ) (H,) S,, fe) dt. 

We prove in section^below that this term has a vanishing contribution, provided 
Tq is large enough. This easy result relies on a non-stationnary phase argument 
in time, recalling that Ugi^t) = exp{—itHg/e) and the energy is larger than 
(5 > 0 . 

• The contribution of very large times, close to the zero energy is 

1 / \ 

-j {u,{t)xsiH,)S„f,)dt. 

We prove in section 0 that this term has a vanishing contribution as £ ^ 0. 
To do so, we use results proved by X.P. Wang |Wa| : these essentially assert 
that the operator (a;)“® Ue{t)xs{He) {x)~^ has the natural size (t)“® as time 
goes to infinity, provided the critical zero energy is non-trapping. Roughly, 
the semiclassical operator Ue{t)xs{He) sends rays initially close to the origin, 
at a distance of the order t from the origin, when the energy is non trapping. 
Hence the above scalar product involves both a function Ue{t)xs {He) Se that is 
localized at a distance t from the origin, and a function fg that is localized at 
the origin. This makes the corresponding contribution vanish. 

The most difficult terms are the last two that we describe now. 

• The contribution of large times, close to the zero energy is 

^ £ e-“'* {Ug{t)xs {He) Se, fe) dt. 

The treatment of this term is performed in section El It is similar in spirit to 
(though much harder than) the analysis performed in the previous term: using 
only information on the localization properties of Ue{t)xs {Hg) Sg and 4>g, we 
prove that this term has a vanishing contribution, provided Ti is large enough. 
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To do so, we use ideas of Bouzouina and Robert |BR| . to establish a version 
of the Egorov theorem that holds true for polynomially large times in e. We 
deduce that for any time Ti < t < the term Ue{t)xs {He) Se is localized 
close to the value at time t of a trajectory shot from the origin. The non-trapping 
assumption then says that for Ti large enough, Ue{t)xs (He) Se is localized away 
from the origin. This makes the scalar product {Ue{t)xs {He) Se,4>e) vanish 
asymptotically. 

• The contribution of moderate times close to the zero energy is 

\ (1 - x) e-“=‘ (Ue{t)xs {He) Se, ^e) dt. 

This is the most difficult term: contrary to all preceding terms, it cannot be 
analyzed using only geometric informations on the microlocal support of the 
relevant functions. Indeed, keeping in mind that the function Ue{t)xs {He) Se is 
localized on a trajectory initially shot from the origin, whereas 4>e stays at the 
origin, it is clear that for times TqE < t <Ti, the support of Ue{t)x5 {He) Se and 
(pe may intersect, due to trajectories passing several times at the origin. This 
might create a dangerous accumulation of energy at this point. For that reason, 
we need a precise evaluation of the semi-classical propagator Ue{t), for times up 
to the order t ^ Ti. This is done using the elegant wave-packet approach of M. 
Combescure and D. Robert |(Tio| (see also |Ro| . and the nice lecture | r}o2| 'l: 
projecting Se over the standard gaussian wave packets, we can compute Ue{t)Se 
in a quite explicit fashion, with the help of classical quantities like, typically, 
the linearized flow of the Hamiltonian — n^{x). This gives us an integral 
representation with a complex valued phase function. Then, one needs to insert 
a last (small) cutoff parameter in time, denoted 0 > 0. For small times, using 
the above mentioned representation formula, we first prove that the term 

^ (1 - x) {Ue{t)xs {He) Se, pe) dt, 

vanishes asymptotically, provided 6 is small, and Tq is large enough. To do so, 
we use that for small enough 6, the propagator Ue{t) acting on Se resembles 
the free Schrodinger operator exp {it[/S.x/2 -\- n^(0)]). In terms of trajectories, 
on this time scale, we use that Ue{t)Se is localized around a ray that leaves 
the origin at speed n(0). Then, for later times, we prove that the remaining 
contribution 

(Ue{t)XS {He)Se,cPe) dt, 

is small. This uses stationary phase formulae in the spirit of EEH, and this 
is where the transversality assumption (H) pageES enters: trajectories passing 
several times at the origin do not accumulate to much energy at this point. 

We end up this sketch of proof with a figure illustrating the typical trajectory 
(and the associated cutoffs in time) that our analysis has to deal with. 
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2.2 Notations used in the proof 

Throughout this article, we will make use of the following notations. 

• Semi-classical quantities 

The semi-classical Hamiltonian and its associated propagator C4(t) have 
already been defined. We also need to use the Weyl quantization. For a symbol 
a{x,^) defined on its Weyl quantization is 

(Opr(a)/) (i) := ^ 

Throughout the paper, we use the standard semi-classical symbolic calculus, 
and refer, e.g., to |DSj or |Ma| . In particular, for a weight m{x,^), we use 
symbols a{x,^) in the class S(rn), i.e. symbols such that for any multi-index a, 
there exists a constant Ca so that 

|5“a(x,e)| < Cam{x,0, \/{x,0 £ 

The notation a ^ ^ £^ak means that for any N and any a, there exists a 
constant CN,a such that 

/ N 

a(x,C) -^e^akix,^) 

\ k=Q 

• Classical quantitities 

Associated with the Hamiltonian JI(x, — n'^{x), we denote the Hamil¬ 

tonian flow 

= {Xit,x ,^), E{t,x,^)), 

defined as the solution of the Hamilton equations 

^X{t,x,^) = E{t,x,0, X{0,x,O=x, 

^E{t,x,0 = iX{t,x,^)), E{0,x,O=i- (2-8) 


< CN,a £^~^^rn{x,^), V(x,^) € 
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(2.9) 


These may be written shortly 

s'®**-*'«> = •' 

where J is the standard symplectic matrix 


J = 


0 Id 


-Id 0 

The linearized flow of d) is denoted by 

DHt,x,0 

It may be decomposed into 

A{t,x,0 B{t,x,i) 


' C(t,x,C) D{t,x,0 J ’ 

where the matrices A{t), B{t), C{t), and D{t) are, by definition 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


A{t,x,^) 

C{t,x,0 


DX{t,x,^) DX{t,x,0 

Dx ' 

D^i,t,X,0 DE{t,x,0 

Dx ’ ■ 


Upon linearizing the matrices A{t), B{t), C{t), and D{t) clearly satisfy 

the differential system 


^A{t,x,^) = C{t,x,^), 

3 

—C(t,x,0 = ^^iX(t,x,0) A{t,x,^), 
together with 

^B(t,x,0 = D{t,x,i), 

3 

—D{t,x,i) = -^{X{t,x,i)) B{t,x,i), 


A(0,a;,^) = Id, 

C(0,x,C) = 0, 

B(0,a:,C) = 0, 


£)(0,x,C) = Id. 


(2.13) 


(2.14) 


In short, one may write as well 

^^F{t,x,i) = J F{t,x,0- (2.15) 

A last remark is in order. Indeed, it is a standard fact to observe that the matrix 
F(t,x,^) is a symplectic matrix, in that 

Fit,x,0'^JF{t,x,0 = J. (2.16) 

for any (t, x, ^). Here, the exponent T denotes transposition. Decomposing F{t) 
as in ifm . this gives the relations 

A{tfC{t) = C{tfA{t), B{tfD{t) = D{tfB{t), 

A{tfD{t) - C{tfB{t) = Id. (2.17) 
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These can be put in the following useful form 

{A{t) + {C{t) + iD{t)) = {C{t) + iD{t)f {A{t) + iB{t)) 

{C{t) + iD{t)f {A{t)-iB{t)) (2.18) 

- {A{t) + iB{t)f {Cit) - iD{t)) = 2ilA. 

These relations will be used in section 0 

3 Small time contribution: the case 0 <t <Tq£ 

In this section, we prove the following 

Proposition 1. We use the notations of section\^ The refraction index is 
assumed bounded and continuous. The data S and (j) are supposed to belong to 
S(MA). Then, the following holds: 

(i) for any fixed value of Tq , we have the asymptotics 

e/ ^ (^) dt 

^ + dt. (3.1) 

(a) Besides, there exists a universal constant Cd depending only on the dimen¬ 
sion, such that the right-hand-side of m satisfies 

i / ^ (^) (*^(^x/2 + ^^(O))) S, (j)) dt - 

<CdTo-^/‘^+^ —> 0. (3.2) 

To —>oo 

Proof of proposition 0 

Part (i) 

In order to recover the limiting value announced in EU, we first perform the 
inverse scaling that leads from to u® (see EU). We rescale time t by a factor 
£ as well. This gives 

= j X {Ue{et)Se,(j)e) dt 

= I X dt. 

We now let 

w'^(t, x) := exp (it (A/2 + n^{ex))) S{x). 

The function w®(t, a;) is bounded in L°° (R; (R‘^)), and it satisfies in the 

distribution sense 

idiw^it, x) = —-Aa;W®(t, x) — n^(ex)w®, w^(0, x) = S{x). 
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These informations are enough to deduce that there exists a function w(t,x) € 
L°° (R; such that a subsequence of x) goes, as e ^ 0, to w(t, x) 

in L°° (R; (R^^)) - weak*. On the more, the limit w(t,x) obviously satisfies 

in the distribution sense 

idt'w{t, x) = x) — n^(0)w, w(0, x) = S{x). 

In other words 

w(t) = exp {it (A/2 + n^(0))) S{x). 

Hence, by uniqueness of the limit, the whole sequence w’^{t,x) goes to •w{t,x) 
in L°° (R; (R'^))-weak*. This proves (ITIII and part (i) of the proposition. 

Part (a) 

This part is easy and relies on the standard dispersive properties of the free 
Schrodinger equation. Indeed, we have 


(exp {it (Aa;/2 + n^(0))) S, (/) 


< 


exp{it{Ax/2 + n'^{0))) S ||(/i||z,i 


\\sh^ II^IUi, 


(recall that S and (j) are assumed smooth enough to have finite norm), for 
some constant Cd > 0 that only depends upon the dimension d. This, together 
with the integrability of the function at infinity when d > 3, ends the 

proof of iTT^ . □ 


4 Contribution of moderate and large times, a- 
way from the zero energy 

In this section we prove the (easy) 

Proposition 2. We use the notations of section n The index is assumed 
to have the symbolic behaviour n~m . The data S and (p are supposed to belong 
to L^(R‘^). Then, there exists a constant Cs > 0, which depends on the cutoff 
parameter S, such that for any e < 1, and Tq > 1, we have 

] (1 - X) ( (1 - X5 {He)) US)Se, 4) dt 


Proof of proposition 121 

The proof relies on a simple non-stationary phase argument. Indeed, this term 
has the value 

^ (1 - X) ( (1 - X5 {He)) exp S^, pe) dt. 
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Hence, making the natural integrations by parts in time, we recover the value 


/{l-XS (H,)) 




A direct inspection shows that this is bounded by 


^ + oo 


Ce^ s-^ l!5|U. ||<^IU. 


<Ce^ S 


2 r-3 


|W3,< 


53 

W 

1 


(1-x) I I e 


Tn ^£2 Tc 


Tn £ 


+ a? + a? 


oe 


dt 


□ 


5 Contribution of large times, close to the zero 
energy: the case t > 


In this section we prove the following 

Proposition 3. We use the notations of section\^ The index is assumed 
to have the symbolic behaviour TTM) . The Hamiltonian flow associated with 
Sfl/2 — n^(x) is assumed non-trapping at the zero energy level. Finally, the data 
S and (p are supposed to belong to 5(R'^). Then, for any <5 > 0 small enough, 
and for any k > 0, there exists a constant Ck ,5 depending on k and 6, so that 


1 

e 



e-“=* (Uflt)xs {He)S„Pe) 


dt 


^ Ck^S £■ 


(5.1) 


The proof relies on the dispersive properties of the semi-classical propagator 
Uflt), inherited from the ones of the classical flow d)(t). More quantitatively, we 
use in this section a Theorem by X.P. Wang |W^ . that we now state. Our index 
of refraction n^{x) is such that n^{x) lies in C'°°(K'^), and it has the symbolic 
behaviour 

n2(x) = - V{x), with |a“P(a:)| < 

(the case 0 < p < 1 is the long-range case, and the case p > 1 is the short-range 
case, in the terminology of quantum scattering). On the more, the trajectories 
of the classical flow at the zero energy (i.e. on the set {(x, f) G s.t. /2 — 
n^{x) = 0}) are assumed non-trapped. It is known inni that this non-trapping 
behaviour is actually an open property, in that 

there exists a i5o > 0 such that for any energy E 

satisfying |if| < i5o, the trajectories of the classical flow (5.2) 

at the energy E are non-trapping as well. 
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Under these circumstances, it has been proved in |Wa| that for any real s > 0, 
and for any ?7 > 0, the following weighted estimate holds true, 

VtGR, (5.3) 

provided the cutoff in energy S satisfies (5 < (5o, i.e. provided we are only looking 
at trajectories having a non-trapping energy. This inequality holds for any test 
function /, and for some constant Cs^n,s depending only on 6, r] and s. In the 
short-range case (p > 1), one may even take p = 0 in the above estimate. Note 
that |Wa| actually proves more: in some sense, the non-trapping behaviour of 
the classical flow is equivalent to the time decay (ESI). We refer to the original 
article for details. We are now ready to give the 

Proof of propositional 

Taking S < Sq, we estimate, using 115.311 . 


/ -I-OO 


p + OO 


< 


{x)-^Ue{t)xsiHe)S4L^ ll(x)>e||L= dt 


dt 


1 r+°^ n. 

<Cs,,,s ||(x)*5,(a:)|U2 ||(a:)«4|U2 

Hence, taking s large enough, and ij small enough, e.g. s = 2 -|- 2 /k, ry = 1, we 
obtain an upper bound of the size 


£ ||(a;)®S'(a:)||L2 ||(a;)®0||L2. 

Here we used the easy fact that ||(a:)®/e{a;)||i 2 < ||(a;)®/(x)||i 2 , when £ < 1, 
together with |j (a;)®5'(a;)||i2 < oo, and similarly for cj). □ 


6 Contribution of large times, close to the zero 
energy: the case Ti < t < 

To complete the analysis of the contribution of “large times” and “small ener¬ 
gies” in (EU that we begun in section 0 there remains to estimate the term 

I £ (1 - X) (J£j (X5 {H,) 4^,) dt. (6.1) 

In this section, we prove the 

Proposition 4. We use the notations of section H The index is assumed 
to have the symbolic behaviour Trm with > 0"^. The Hamiltonian flow 
associated with ^ j^ — vfl'ix) is assumed non-trapping at the zero energy. Finally, 
the data S and 4> are supposed to belong to 5(R'^). Then, for 5 > 0 small enough, 

'^The assumption is crucial, see Lemmabelow. It ensures that the wave Ue{t)Se 
propagates with a uniformly non-zero speed, at infinity in time t. 
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there exists a Ti ((5) depending on S such that for any Ti > Ti (5), we have for 
K small enough, 

\jl (1 - (^) 

< Ck ,5 e, as £ ^ 0, (6.2) 

for some constant Ck,s that depends upon k and 6. 


The idea of proof is the following: the functions Sg and (j)g are microlocally 
supported close to points {xo,£,o) & such that xq = 0 (due to the concen¬ 
tration of both functions close to the origin as e ^ 0). On the more, using 
the Egorov Theorem, one may think of the time evolved function Ug(t)Sg as 
being microlocally supported close to points (X(t; xq,■ fo),2(t; ccq,^ o)) that are 
trajectories of the classical flow, with initial data {xo,fo) such that xq = 0. 
Using the non-trapping assumption on the classical flow, we see that for large 
times t > Ti with Ti large enough, the trajectory X(t;xo,^o) with xq = 0 is 
far away from the origin. Hence the microlocal support of Ue{t)Ss and (fe do 
not intersect, and the factor EU should be arbitrary small in e as e ^ 0. 

The difficulty in making this last statement rigorous lies in the fact that we 
need to use the Egorov Theorem up to (polynomially) large times of the order 
t ~ This difficulty is solved in Lemma 0 below. Indeed, upon adapting a 
recent result of Bouzouina and Robert |BR| we give remainder estimates in the 
Egorov Theorem that hold up to polynomially large times (logarithmic times 
are obtained in the context of [BB j 1. This is enough to conclude. 


6.1 Proof of proposition m 

The proof is given in several steps. 

First step.' Preliminary reduction 

In this step we quantify the fact that the functions involved in the scalar 
product in are microlocalized close to the zero energy = n^(x) (in 

frequency) and close to the origin x = 0 (in space). To do so, we simply write, 
using the fact that S and (f belong to 5(M‘^), 

= xsi\x\)(l)eix) + Osi£°^) in L^(K‘^), 

and similarly for S^. This means that for any integer N, there exists a Cn ,5 > 0 
that depends on N and 6, such that Wfieix) — X 5 (|x|)(/)e(x)||j;^ 2 (R<i) < Cn£^. As 
a consequence, we may rewrite the contribution we are interested in as 

^ (1 - x) (X5(|2;|) X5 (He) Ue(t)x5i\x\) S'e, (jje) dt 

up to an Os{£°°). There remains to bound the above term by 


< Il'S’ellia II'/'eIIl^ X - / 

£ Jt, 


xs (|a;|)X5 {He)Ue{t)xs (kl) 


C{L'^ 


£ Jti 


xs {\x\)xs {H^)Ue{t)xs (kl) 




dt , 


dt 


(6.3) 
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up to an Os{e°°). Our strategy is to now evaluate the operator norm under the 
integral sign. This task is performed in the next two steps. 


Second step: symbolic calculus 

In view of (lOl . our analysis boils down to computing, for any Ti < t < e 
the operator norm 


xs (kl)X5 {He) Ue{t)xs (|a;|) 


2 

C{L^)' 


Expanding the square, this norm has the value 


XS (|a::|) U*{t)xs (He) Xs (|a:|) {He) Ue{t)xs (|a:|) 


C{L^)' 


(6.4) 


Now, and for later convenience, we rewrite the above localizations in energy and 
space, as microlocalisations in position and frequency. 

Using the functional calculus for pseudodifferential operators of Helffer and 
Robert |HR| (see also the lecture notes m and [Map , there exists a symbol 
Xs{x, such that 


(He) = Op^iXs) + 0(£°“) in £(L2). 


The symbol Xs{x,^) is given by a formal expansion 

( 6 . 5 ) 

fc >0 

where the expansion (16.51) holds in the class of symbols that are bounded to¬ 
gether with all their derivatives. Furthermore, the principal symbol of Xs is 
computed through the natural equality 

Xs^°\x,C) = ■ 

Finally, the explicit formulae in [DSj give at any order A: > 0 the following 
information on the support of the symbols Xs^^\ 

supp Xs^'"'^ C {|C^/2 - n^{x)\ < <5}. 


Hence (16.411 becomes, using standard symbolic calculus, 


X5(kl) U*{t) [Op“(A’5(a;,^) ‘^x}{\A) tt0)] Ueit)xs{\x\) 


C(L^) 


,( 6 . 6 ) 


up to an Os{e°°) (Here we used the uniform bound \\Ue{t)\\c{L‘^) < !)• Let 
us define for convenience the following short-hand notation for the symbol in 
brackets in H6.6II : 


hs{x,C) '■= Xs{x,0 HXs (kl) {iXsix,^- 


The only information we need in the sequel is that bs admits an asymptotic 
expansion bs = where each bg^'^ has support 

supp C {|a:| < (5} (b {|^^/2 — n^(x)| < (5} =: F(5). 
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This serves as a definition of the (compact) set E{6) in phase space. In the 
sequel, we summarize these informations in the following abuse of notation 

svcpp bs C E{S). (6.7) 

The remainder part of our analysis is devoted to estimating 

X5(|a;|) U*(t)Op'^ibsix,0)Ueit)xsi\x\) , 

C{L^) 

and the hard part of the proof lies in establishing an “Egorov theorem for large 
times”, to compute the conjugation U*{t)Op'^{bs{x,^))Ue{t) in 16.411 . 


Third step: an Egorov theorem valid for large times - End of the proof 
Now we claim the following 


Lemma 5. We assume that the refraction index has the symbolic behaviour 
xrm with > 0^. We also assume that the zero energy is non-trapping for 
the flow. Take the cutoff parameter in energy 5 small enough. Then, 

(i) Let $(t, X, f) be the elassical flow associated with the Hamiltonian ^ jl — 
n^ix). Let F(t,x,f) he the linearized flow. For any multi-index a, and for any 
(small) parameter r] > 0, there exists a constant Cs^\a\,rj such that for any initial 
datum {x,f) € E{5) = {|x| < d} H {|^^/2 — n^(x)| < d}, we have 


yt e R, 


d^F{t,x,0 

dix,f)‘^ 


<CsM,v (t)(i+’’)(i+l“l)+2|“l. 


( 6 . 8 ) 


In other words, the linearized flow has at most polynomial growth with time, 
(a) As a consequence, for any time t, there exists a time-dependent symbol 

hs(t,x,^) - ^e%s^^\t,x,^), 

k>0 


such that the following holds: there exists a number cs > 0 such that for any 
N > 0, there exists a constant such that 


N 


u:{t)0pf {bs) US) - Op“ ^ eSs 


(k) 


\k=0 


(6.9) 


C(L^) 


Again, the error grows polynomially with time, and we have some control on the 
dependence of the estimates with the truncation parameter N. 

(Hi) Moreover, we have the natural formulae 
hs''°\t,x,f,) = bs S{t,x ,^)), 

and, for any k > 0 we have the information on the support 

supp a;,^) C {(x,^) G R^^^ s.t. ^{t,x,f,) e E{S)}. 


^The assumption > 0 is crucial, see 
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We postpone the proof of Lemma 13 to paragraph |^1 below. We first draw its 
consequences in our perspective. 

Leaving TV as a free parameter for the moment, we obtain 


£(L2) 


£(L2) 


X<5 (|a:^l) U*{t)Op^{bs{x,0)US)xs (|a;|) 

X5(kl)Op^ I j X5(|a;|) 

\fc=o / 

+ Os 

Op“ (x^(N) tt ) ttX6(N) 


\k^0 


C{L^) 


Os . 


Now, part (iii) of Lemma [S] and standard symbolic calculus indicate that the 
above symbol has support® in 


fe^o (|a:^l) nsupp 

C {(a:,^) s.t. |a;| < 6, and G if((y)}. 


The non-trapping condition (and more precisely estimate Iti.lOII below) allows 
in turn to deduce that this set is void for t large enough. Hence, up to taking a 
large value of Ti, Ti > Ti (5) for some Ti (^), we eventually obtain in Ifi.lfll . 


< 


'Ti 


£(L2) 


dt 


(|a;|)X5 {H^)Ue{t)X8{\x\) 

J^Os dt < Os (e(A^-i)/2-c..ivV2^ < 


for K small enough (and TV = 4 will do). This ends the proof of propositional 


6.2 Proof of Lemma m an Egorov theorem for polynomi- 
ally large times 

In view of the above proof, we are left with the task of proving the large time 
Egorov theorem of LemmaEl To do so, we follow here closely ideas developped 
in |BR| in a slightly different context. Part (iii) of the Lemma is proved in [BR| . 
so we will skip this aspect. The implication (i) =J> (ii) in Lemma 0 which we 
prove below for completeness, is also essentially proved in [BR] . Our main task 
in the sequel turns out to be the proof part (i) of the Lemma. 

The proof is given in several steps. 

First step: estimates on the flow <i>(t,x,^) 

®we make here the same abuse of notation than in lEJ. 
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In this step, we prove that for small enough a (5, there is a time T{S), depending 
on 6, such that for any initial datum (cc,^) of phase-space in the set E{5) = 
{|a;| < 5} n {|^^/2 — n'^{x)\ < 5} ('see l6.7ll '). one has 

Vt>T((5), \X{t,x,0\>Cst, (6.10) 

for some constant Cs > 0 that depends on 6, that is however independent of both 
time t and the initial point under consideration. The proof is standard 

and uses the information > 0. 

First, the non-trapping condition implies that for any large number R' > 0, 
and for any initial point (x,^) S E{d), there exists a time T{R',x,^) such that 

yt>T{R',x,0, \X(t,x,^)\> R'. 

By continuous dependence of the flow X(t,x,^) with respect to the initial data 
(a;,^), and compactness of the set E{6), there is a time T{R',6), that now 
depends upon R' and 6 only, such that for any initial point (a;,^) € E(S), there 
holds 

Vt>T(R'), lX(t,x,OI > Ji'- 

In other words, the trajectory X(t, x, goes to infinity as time goes to infinity, 
uniformly with respect to the initial datum (a;,^) € E(S). 

Second, we get estimates for the standard “escape function” of quantum and 
classical scattering, namely the function X{t) ■ S(t). We compute 

^ iX{t,x,0 ■ E{t,x,^)) = 2 _ ^2 

+ 2n'^ {X{t,x,^)) +X{t,x,^) ■ {X{t,x,^)) 

= 2 (y - n^ix)^ + 2n? (X(t, x, ^)) -f X(t, a;, • Xv? {X{t,x,0) 

(thanks to the conservation of energy) 

t^2(^-n^(x))+2n^, 

uniformly with respect to the initial datum (x,^) G E{S). Hence, using the fact 
that > 0, and taking a possibly smaller value of the cutoff parameter 5, we 
obtain the existence of a constant Cs > 0, and another time T{S), such that 

Vt > T((5), X{t,x,^) ■E{t,x,^) > Cs t. (6.11) 

Using the fact that ^ X‘^{t,x,^)) = X(t,x,^) ■ S(t,x,^), we deduce the 

desired lower bound 

Vt>T(J), l{x\t,x,0-X^{T{S),x,0)>Cs^^. 


Second step: estimates on the linearized flow F{t^x,^). 

One first proves the estimate (lira in the case a = fl = 0. By its very definition 
(ICTTl . the linearized flow 


E{t,x,^) 


A{t,x,0 

C{flx,i) 
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B{t,x,^) 

D{t,x,0 



satisfies (see (ITTHl . (EHU) the differential system 


^(0,x,C) = Id, 
C(0,a:,e) = 0, 


^C{t,x,0 = D'^n^ {X{t,x,^)) A{t,x,^), 
together with 

^B(t,x,0 = D(t,x,0, 

X, C) = -D^n^ iX{t, X, 0) B(t, x, C), 


( 6 . 12 ) 


B{0,x,O = 0 , 

11(0, a;, C) = Id. (6.13) 


Here, the notation D'^n?{x) refers to the Hessian of the function n'^{x) in the 
variable x. Due to the assumption (tmii on the behaviour of n^(x) at infinity, 
we readily have 

\D^n‘^{x)\ < C {x)-P-\ 

for some constant C > 0, independent of x. This, together with the previous 
bound (lOl on the behaviour of the flow X{t,x,^) at infinity in time, gives 
the estimate 


\D^n^{X{t,x,0)\<Co (6.14) 

for some constant Co > 0 which is independent of time t > 0, and of the point 
(x,^) in phase-space. We are thus in position to estimate A{t) and C{t) using 
(IHdH) . Integrating in time, and setting 

£(t) := |DV(X(t,x,e))| (6.15) 

for convenience, we obtain (dropping the dependence on (x,^) of the various 
functions), 

(t — s) e(s) |H(s) — Id| ds-I- f {t — s)e(s)ds, (6.16) 

Jo 

|C(t)| < f e(s) |A(s)| ds. (6-17) 

Jo 

Choose now a constant C*, and define the time t* as 


\m - idi < 


C := sup{t > 0 s.t. \A{t) — Id| < C*(t)^^''}. 


We prove that C = +oo, provided C* is large enough. Indeed, for any time 
t using 16.1611 together with the decay 16.1411 . we have 


\A{t) - Id| < CoC* 



-s)(s)-'’-^+''ds<CoC*t 



-P-^+^ ds 


< C^C-tCrf t 

(for some constant C^ > 0, provided rj > 0 satisfies -q < p 12) 
<C. 

(provided t is large enough, t > T{Co, Cn), for some T{Co, Crj) 
that only depends on Cq and Crj). 
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On the other hand, we certainly have |A(t) — Id| < for bounded values 

of time t < T{Co, Crf), provided C* is large enough. Hence O = +oo. Inserting 
this upper-bound for A in (lOTIl gives 

for some > 0, provided 77 > 0 is small enough. We may estimate B{t) and 
D[t) in the similar way. The analysis is the same, and starts with the formulae 

\B{t)\ <t+ f {t- s) e(s) |H(s)| ds, 

Jo 

\D{t)\ <1+ [ e{s) |i?(s)| ds. 

Jo 

We skip the details. At this level, we have obtained the bound 

\Fit,x,0\<C^ 

for any (small enough) ij > 0, and a constant Crj independent of 


Third step: estimates on the derivatives of the linearized flow 
Let now a be any multi-index. We prove (Oil by induction on |a|. 
any p > 1 


Mp(t) 


sup sup 
|/3|—n 


d(x,iY 


Define, for 


We have proved in the second step above that 




Assume that for some integer pq, the estimate 

Mp{t) < (i)P(i+r/)+2(p-i)^ 


has been proved for any p < po- We wish to prove the analogous estimate 
for Mpg^i. Take any multi-index a of length |a| = po- From now on, we 
systematically omit the dependence of the various functions and derivatives 
with respect to (x,^), and write d°‘F{t), d^Fd instead of 5“F(t,x,^)/5(x,^)“, 
d°‘H{x,^)/d{x,^)°^ and so on. Upon differentiating a times the linearized equa¬ 
tion 112.1511 on F, we obtain, 

dt id°‘F{t)) = J XI ( 5 ) . (6.18) 

l3<a \ ^ 2 

In order to make estimates in dnuHi), we first need to write the Faa de Bruno 
formula as 

{D^H o $(t)) = ;3! ^ (d-D^H) o cl>(t) x H ■ 


Here [3 G 7 G and ( G are multiindices, and m associates to each 
multi-index ( G N^'^, another multi-index m(() G N^'^. Also, the above sum 
carries over all values of 7, m, and ( such that 

Xto(c) = 7 . 

C C 
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Finally, when |/3| > 1, the above sums carries over 7 ’s and C’s such that jyj > 1 
and Id > 1- All this gives in (I6.18II . 


dt [d^F{t)) = J ^ /3! 

0<(y. 


■' 7,m 





Hence, putting apart the contribution stemming from /3 = 0, we recover 


dt {d^F{t)) = J D^H ($(t)) (d^Fit)) + (6.20) 


where the remainder term Ra{t) is estimated by 


\Ra{t)\ 

0#/3<a7,™ C 

0y^l3<a‘y,m (, 

for some constant C\a\ > 0 that depends on |a|. The last line uses the fact 
that supj, for some constant C^. Using the inductive 

assumption, we recover 

07 ^/ 3 <q; 7 ,ni 

X i™(oi 

C 

<C\a\n E 

0^^<a 

— (7|^| ^ ^ ^^^(l+'/)(l+|a-/5| + l/3|)+2(|a-/3| + |/3|-|7|) 

0^/3<q 

< qal.r, (i)(l+r))(l + |a|)+2(|a|-l) _ 

Here we used the constraints together with the information jyj > 1. Using 

Lemma El below in equation satisfied by d°‘F, we obtain. 


Hence 

This ends the recursion. 

Fourth step: A Gronwall Lemma for solutions to the linearized Hamilton equa¬ 
tion 

The preceding step uses the following 
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Lemma 6. Assume the function G{t,x,ff) satisfies the differential equation 

=J-D^H mf x,f))- Git, x,0+O {{t)>') , (6.21) 

G(0,a:,C) = 0, 

where the O is uniform in {x,f,). Then, G satisfies the uniform estimate 

G{t,x,C} = O ((t)^+^) . 

Proof of Lemma El 

Decompose G(t) = G{t,x,f,) as 


G{t) = 


Ac it) Bait) \ 
Gait) Doit) )■ 


Then, equation for G writes 


^^AGit) = GGit) + Oiit)^), 

^CGit) = dW iXit)) AGit) + O iit)^) , 
together with 

^^BGit)=DGit)+0{it)^), 

^DGit) = iXit)) BGit) + O (it)^) , 
Equations give rise to the estimates 


^g(O) = 0 , 

Gg(0) = 0, 


BGiO) = 0, 


Gg( 0 ) = 0 . 


( 6 . 22 ) 


(6.23) 


I^G(t)|<G 
\CGit)\ < G 



- s) ( e(s) |^g(s)| + (s)^) ds, 



|^g(s)| ds, 


(6.24) 

(6.25) 


where the function e(s) is defined in 16.1511 above. Using e(s) < Go (s) ^ ^ < 
Grj (s)“’'“^ for any small r] > 0 (see (Em), gives in equation Em. 


\AGit)\ <Gr,t is)-^-^ |Ag(s)| ds + G it)^+^ . (6.26) 

^0 

From this it can be deduced that 


I^g(0I < c (t)^+2. 

(for a given constant G*, define indeed t* = sup{t > 0 s.t. \AGit)\ < G* (t)''+^ 
- one deduces from (IH:^ that t* = +00 provided G* is large enough - see (IHT^ 
and sequel for details). Equation 16.2511 then gives 

|GG(t)| < G,, |^g(s)| ds < G, (t)^+i-T 

Jo 

The estimates for Bg and Dg are the same. This ends the proof of the Lemma. 
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Fifth step; adapting the estimates of m 

We now put together the estimates on the linearized flow obtained before, to 
complete the proof of parts (ii) and (iii) of Lemma |S| 

The construction of the symbols (t, a;, in Lemma El is made in an 
explicit way in [Bli j . Part (iii) of LemmaElfollows. Also, the remainder estimate 
(EHI) is a consequence of the above estimates on the linearized flow F{t,x,f) 
and its derivatives, upon adapting the analysis of |BR| . Let us indeed write the 
rough (but simpler) estimate 

|a“F(t,;r,C)| 

corresponding to the special choice p = 1 in lOll . Then, Theorem 1.2 - formula 
(12) of [BR] . 

hs^°\t,x,£,) = bs ($(t,a;,C)), 

together with the Faa de Bruno formula, give for any multi-index a the estimate 

|a“b5(°)(t,x,c)| <qa| 

From Theorem 1.2 - formula (14) of [BRj . we have for any fc > 1 the explicit 
value 


hs^'"\t,x,C) = 


E 

|a|+£—/c+l 
0<£<fc-l 



X 


o $(t 


s, x^f) ds, 


where r(a) is a harmless coefficient whose explicit value is given in [BR] . This, 
together with the Faa de Bruno formula, implies for any fc > 1, the upper-bound 

|5“b5«(t,a;,q| (tyoikH+k^+i)^ 

for some fixed number cq, independent of a and k. Then, using formulae (51), 
together with (52), (54), (97) and (99) of [BRj gives the estimate (|6.9|l . This 
ends the proof of Lemma 0 □ 


7 Contribution of moderate times, close to the 
zero energy 

After the work performed in sectionsElthroughEl there only remains to estimate 
the most difficult term 

^ (1 - x) e-“'‘ {Ueit)xs {He) Se. (fe) dt. 

This is the key point of the present paper. 

The main result of the present section is the following 

Proposition 7. We use the notations of section\^ The index is assumed to 
have the symbolic behaviour n~w . The zero energy is assumed non-trapping 
for the Hamiltonian q/2 — n'^{x). Finally, we need the tranversality condition 


28 




(H) vaae l,V.5l on the trajectories ^{t,x,^) with initial data satisfying x = 0, 
f'^/2 = n^(0). Then, the following two estimates hold true, 

(i) for any fixed value of the truncation parameters 9, Ti and 6, we have 

(0 dt -^0. 

(ii) for 9 > 0 small enough, there exists a constant Ce > 0 such that for any 
e < 1, we have 

(ifc) * (9 ■“ 

<CeTo-‘^/‘^+^ —> 0 . 

To —>+oo 


The remainder part of this paragraph is devoted to the proof of proposition [3 
In order to shorten the notations, we define 


aw := (1 - W (j^) e-"W 

SO that the proof of proposition [3 boils down to estimating 


(7.1) 


- Xe{t){xs{He)Se,Ue{-t)(j)e) dt. (7.2) 

^ JToe ' ' 

The precise value of the cut-off function Xe(t) in the analysis of 17.211 will be 
essentially irrelevant in the sequel. 

Proof of proposition [3 

The proof is given in several steps. As in section |3 we begin with some 
preliminary reductions, exploiting the informations on the microlocal support 
of the various functions. Then, we use the elegant wave-packet approach of 
Combescure and Robert ISi to compute the semi-classical propagator Ue{t) 
in in a very explicit way - see Theorem|3 below: this gives a representation 
in terms of a Fourier integral operator with complex phase, that is very well 
suited for our asymptotic analysis (see also [(Ti.Rj . or the work by Hagedorn 
and Joye EH, |H2| . jH.T| . or by Robinson |R b| . or even the seminal work by 
Hepp m for similar representations - see also Butler |Bt|i. This eventually 
reduces the analysis to stationary phase arguments that are very much in the 
spirit of ITTr^ . and where the tranversality assumption (H) page I35| turns out 
to play a crucial role. 


First Step: Preliminary reduction, projection over the gaussian wave packets 
As in sectionEl(see H6.3II . (16.511 . (I6.7|l l. we may first build up a symbol ao(x,f) G 

suppoo C {|a:| < (5} n {|C/2 — n‘^ix)\ < 5}, (7.3) 


and 

(H,) S„ U,i-t)cf,) = (Op“ iaoix,f)) S„ U,i-t)j,,) + Osie°^). 
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With the notation (lti.511 . we actually have the value ao(a;,^) = Xs{x,^)^xs (|a;|)- 
Therefore, the asymptotic analysis of reduces to that of the expression 


- Mt){Op^{ao)Se,U,{-t)cl,,)dt. (7.4) 

^ JTos '' ' 

Now, to be able to use the wave-packet approach of we need to decompose 

the above scalar product on the basis of the Gaussian wave packets 

:= exp (^-p ■(*“!)) (~ 2^ ) ' 

Each function is microlocally supported near the point {q,p) in phase-space. 
Using the well-known orthogonality properties of these states, i.e. 


{u,v) = (27r£) ^ / dqdp {u,Lp\ p) 

for any u{x) and v{x) in the space and forgetting the normalizing factors 

like TT, etc., we obtain in mil 


1 p 

/ / dtdqdpXs(t)(Op^{ao)Se,iplp) (plp,Ue{-t)(j)e) 

(Ue{t)plp,(l)ey (7.5) 

Before going further, and in order to prepare for the use of the stationary phase 
theorem below, we make the simple observation that the integral dqdp over 
in csi may be carried over the compact set {|a;| < 26} n {|C^/2 — n'^{x)\ < 26}, 
up to a negligible error Os{s°°). For that purpose, take a function Xo{q,p) S 
C'“(M^‘^) such that 


suppxo(9,p) C {|x| < 2^} n {|C^/2 - n^{x)\ < 26} 

Xo( 9 ,p) = 1 on {|a;| < 35/2} n {|^^/2 - n^(a;)| < 3(5/2} . (7.6) 


We claim the following estimate holds true: 


/R2‘i 


dqdp (1 - Xo{q,p)) Op™(ao)</5®_p 


L 2 (R<i) 


= Os{e^). 


Indeed, we have the following simple computation: 

|2 


(7.7) 


Op/'(ao)<Pg,p||^2(Rd) = {OPeMao)iplp,iplp) 

= [ dxd^{ao'iao){x,^) W{iplp){x,^) 

(where W{‘p}^ p) denotes the Wigner transform of <p\^p) 


= £ 


dxdS, (aojlao)(a;,0 exp ( - 


\q-x\^ -b \p-i'\‘ 


and the last line uses the fact that the Wigner transform of ip^^p is a Gaussian. 
Now, using supp(ootlao) C {|a;| < 5} n {|^^/2 — n^(x)| < 6}, together with 17.611 . 
establishes (TT^ . 
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Using this estimate G3, and replacing back the factor Op™(ao) by the 
identity in we arrive at the conclusion 



To E 


Our strategy is to now pass to the limit in the term 



In order to do so, we need to compute the time evolved gaussian wave packet 
Ue(t)ipg p in an accurate way. 

Second Step: Computation ofUe(t)ip^q p - reducing the problem to a stationary 
phase formula 

The following theorem is proved in (see also [E3, EoS) 

Theorem 8. (|CRgJ, lEo] ) We use the notations of section ^ Under as¬ 
sumption irm on the refraction index n^{x), there exists a family of functions 
{pfej (t, g,p, x)}(fcj)gN 2 , that are polynomials of degree at most k in the variable 
X G with coefficients depending on t, q, p, such that for any e < 1, the 
following estimate holds true: for any given value of Ti, and any given integer 
N, we have, for any time t G [0, Ti ], 



where 


QN{t,q,p,x) := 1 + 



/at ;= {1 < j < 2iV - 1, l<k-2j <2N -I, k> 3j} . 


Here, the following quantities are defined: 
• Ag is the dilation operator 



(7.10) 


• %(qt,pt) is the translation (in phase-space) operator 



(7.11) 


• {Qt,Pt) denotes the trajectory 

{qt,Pt) ■■= {X{t,q,p),E{t,q,p)) 


(7.12) 
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• S(t,q,p) denotes quantity 


(7.13) 


S{t,q,p)=J^ (^^+n^{qs)^ rfs _ 

• Ai{F{t,q,p)) is the metaplectic operator assoeiated with the sympleetic matrix 
F(t,q,p). It acts on the gaussian as 

M{F{t,q,p)) (^exp = 

det(A(t, q,p) + iB{t, q,p))-^/‘^ exp . ( 7 . 14 ) 

Here, the square root det{A{t, q,p) + iB(t,q,p))c ' is defined by continuously 
(hence the index c) following the argument of the complex number det{A{t, q,p) + 
iB{t, q,p)) starting from its value 1 at time t = d. Also, the complex matrix 
T(t,q,p) is defined as 

r(t,g,p) = {C(t,q,p) +iD{t,q,p)) {A(t,q,p) + iB{t,q,p))~\ (7.15) 

Remark 

If the refraction index n^{x) is quadratic in x, then formula lEnj is exact, and 
the whole family {pkj} vanishes. This is essentially a consequence of the Mehler 
formula. We refer to for a very complete discussion about the propagators 
of pseudo-differential operators with quadratic symbols. 

In the case when n^(x) is a general function, the polynomials pkj are ob¬ 
tained in EESI using perturbative expansions “around the quadratic case”. We 
refer to m for a very clear and elegant derivation of these polynomials. Let us 
quote that similar formulae are derived and used in m- The idea of consider¬ 
ing such perturbations “around the quadratic case” traces back to |He| , see also 

[HT1 . EH, 1^ - 

The fact that the matrix A{t) + iB{t) is invertible, and r(t) is well defined, 
is proved in [^, see also |R,o2| . It is a consequence of the symplecticity of F{t) 
(see the relations We refer to the sequel for an explicit use of these 

important relations. □ 

In the next lines, we apply the above theorem, and transform formula (EH) 
accordingly. 

On the one hand, we use the Parseval formula in (EH to compute the two 
scalar products. Forgetting the normalizing factors like tt, etc., it gives, e.g. for 
the first scalar product, 

{Se, Vq,p) = f dxdf, exp(ia; • ^/e) S{£,) Pq^p{x) 

= [ dxdf, exp(ia: • (,/e) xi{x)S{() ^(x) -f 0(e°°), 

for any truncation function xi being = 1 close to the origin. On the other 
hand, we use formula (17.911 to compute Ue{t)(pg p in (17.811 . using the short-hand 
notation 

PN{t,q,p,x) := det{A(t,q,p) + iB(t,q,p))/:^/'^ QN{t,q,p,x). 
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These two tasks being done, we eventually obtain in upon computing 

the relevant phase factors explicitly, 


1 

e 


[ Xe{t) (xs {He) Se, Ue{-t)(j)e) dt = Ot^ ,5 (£°°) + 
JTos '' ' 


1 




£.(5d+2)/2 / 

-'To { 


dtdqdpd^dijdxdy Xe{t) exp ( -^{x,y,^,y,q,p,t) 


S{0^*{v)xo{q,p)xi{x,y)PN ( t,q,p. 


y-qt 


(7.16) 


where xi G C'“ is = 1 close to (0,0). Here, the crucial (complex) phase factor 
has the value 


'i’{x,y,^,r],q,P,t) = J ^y+n^(gs)^ ds - p ■ {x - q) + pt ■ {y - qt) 


+ x- ^ — y-rj + i 


.{x-q)'^ T{t){y-qt) ■ (y-qt) 


(7.17) 


Our goal is now to apply the stationary phase formula to estimate (TtTtIi . Ob¬ 
viously, the cutoff in time away from t = 0 in prevents one to use directly 

the stationary phase formula close to t = 0. This is the reason why times close 
to 0 are treated apart in the sequel (see steps four and five below - see also the 
outline of proof given in section P)l. 


Third Step: computing the first and second order derivatives of the phase 
First, it is an easy exercice, using the symplecticity relations (ITTTIl . to prove 
that the matrix T{t) is symmetric and it has positive imaginary part. The 
relation 

Im {T{t){y-qt) ■ {y-qt)) = \ {A{t) + iB{t))~^ {y - qt)\'^, 
implies indeed 

Im^* = |x - qp -b I {A{t) + iB{t))~^ {y - qt)\^. 

Hence we recover the equivalence 


Im = 0 iff y = qt and x = q. 

Second, using the differential system 
A{t), B{t), C{t), and D{t), we prove 

Vq.p (^y +«^(9s)^ ds^ = 

This gives the value of the gradient of 'k 


(7.18) 

satished by the matrices 
A{t)'^Pt - P \ 

B{t)^Pt J ■ 


V 




'^{x,y,^,V,q,p,t) = 


( -p + i + i{x-q) \ 

Pt-y + T{t){y-qt) 

X 


c{t)^{y - qt) + i{q -x) + A{t)'^r{t){qt - y) 

-{x -q) + D{t)'^{y - qt) + B{t)'^T{t){qt - y) 

„2 

\ -y + n'^{qt) + yn'^{qt) ■ {y - qt) + Pt ■ ^{t){qt - y) J 


(7.19) 
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This computation is done up to irrelevant O ((y — qtY' + {x — terms. 
These observations allow to compute the stationary set, defined as 

M := {(x, y, 77 , q,p, t) S IR“x]0, +oo[ 

s.t. Imd7 = 0 and Vx,y,i,ri,q,p'^ = 0}. (7.20) 

Note (see above) that we exclude the original time t = 0 in the definition of M. 
In view of mu . the set M has the value 

M = {{x,y,^,q) s.t. x = y = q = 0, £,=p} 

n |(p,77,t) s.t. ^ = n^(0), gt = 0, pt = 7?! . (7.21) 

Note that the second set reads also, by definition, 

|(p,77,t) s.t. y = 77,2(0), X(t,0,p) = 0, S(t,0,p) = Tyj . 


Last, there remains to compute the Hessian of d' at the stationary points. A 
simple but tedious computation gives, for any point {x,y, ^,r], q,p,t) G M, the 
value 


D 


2 




{x ,y ,7],q,p,t) ^M 


ild 

0 

Id 

0 

-ild 

-Id 

0 

0 

Tt 

0 

-Id 

Ct-TtAt 

Dt-TtBt 

vA(o) 







-rtv 

Id 

0 

0 

0 

0 

0 

0 

0 

-Id 

0 

0 

0 

0 

0 

-ild 

1 

0 

0 

0 

-CjAt+ild 

-CjBt 

-cjv 





AA'^ TtAt 

AA'^TtBt 

AAjrtr, 

-Id 


0 

0 

Id-DjAt 

-DjBt 

-Djv 





ABjTtAt 

ABjFtBt 

AB^rtV 

0 

Vn^(O)'^ 

0 

0 

-V^Ct 

-v'^Dt 

-AVrt=(0) 


-r7'^rt 



+rPTtAt 

+ Art St 

+r;'^rt77 


\ 




Here we wrote systematically At, Bt, etc. instead of A{t), B{t), etc. The above 
matrix is symmetric, due to the relation The very last computation we 

need is that of Ker at stationary points. The value of clearly shows 

that 

Ker (772^1^) = {(X, r,S,77, Q,P,T) s.t. A = y = Q = 0 , S = P, 

-77 + (A - TtBt)P + T(Vr 72 ( 0 ) - ^77) = 0 , 
{-Cj + AjTt)BtP + T{-Cj + AjTt)y = 0 , 
{-Dj + BjTt)BtP + T{-Dj + BjTt)y = 0 , 

Tf{-Dt + TtBt)P + r77T(-V772(0) + rt77) = 0}. 
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Hence, using Dj — BfVt = {At + iBt) together with Cj — AjVt = —i{At + 
iBt)~^, and {At + iBt)~^’'^ + TtBt = Dt (see 112.1811 h we obtain 

Kev {DH\m) =[{X,Y,^,H,Q,P,T) s.t. X = Y = Q = 0,E = P, 

and 77 '^H = 0, BtP + T ?7 = 0, H = DtP + TVn^{0) = o|. (7.22) 


Fourth Step: Application of the stationary phase Theorem - proof of part (i) 
of proposition^ 

In this step, we formulate the main geometric assumption on the flow $(t, x, ^), 
that allows for the proof that the contribution in rmiii vanishes asymptotically. 

(H) Transversality assumption on the flow 

We suppose that the stationary set 

M = {x = y = q = 0,f = p}n |y = n2(0), X{t,0,p) = 0, E{t,0,p) = 7 ?| 

is a smooth submanifold o/]R®'^x] 0 ,+oo[, satisfying the additional constraint 

k := codimM > 5c? + 2. (7.23) 

We also assume that at each point m = {x,y,f,ri,q,p,t) G M, 
the tangent space of M at m is 

TmM = ^{X,Y,E,H,Q,P,T) s.t. X = Y = Q = 0,E = P, 

and rj^H = 0, BtP + Ty = Q,-H + DtP + TVn'^{Q) = o|. (7.24) 

In other words, we assume that T^M is precisely given hy 
linearizing the equations defining M. 

First Remark 

We show below examples of flows satisfying the above assumption. It is a 
natural, and generic, assumption. Note in particular that the assumption on 
the codimension is natural, in that the equations defining M give (roughly) 
4c? constraints on {x,y,q,f), one constraint on rj, and again 2d constraints on 
the momentum p, the solid angle r]/\r]\, and time t. Hence one has typically 
A: = 6 c? + 1. □ 

Second remark 

Equivalently, the above assumption may be formulated as follows. The set 

M := {{p,V,t) s.t. y = 71^(0), X{t,0,p) = 0, 5(t,0,p) = ry} 

is assumed to be a smooth submanifold of satisfying the additional con¬ 

straint codim A4 > c? -|- 2 , and whose tangent space is given by 

{{P, H, T) s.t. rj^H = 0, BtP -h Try = 0, DtP + TV71^(0) -H = Q}. 

Note in passing that the conservation of energy allows to replace the requirement 
7 y ^/2 = 77 ^( 0 ) by the equivalent p ^/2 = 77 ^( 0 ) in the definition of M. □ 
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Third Remark 

Provided M is a smooth submanifold with tangent space given upon linearizing 
the constraints, its codimension anyhow satisfies 

codim M > 5d + 2. 


Equivalently, provided is a smooth submanifold with the natural tangent 
space, its codimension anyhow satisfies 

codim Af > d + 2. 


As a consequence, the analysis given below (see ( 17.2711 1 establishes that {w^,(j)) 
is uniformly bounded in e. This fact is not known in the literature. □ 

Under assumption (H), we are ready to use the stationary phase Theorem in 
(Trroii . at least for large enough times t (recall that the very point t = 0 is 
excluded from the definition of M above). Indeed, assumption (H) precisely 
asserts the equality 

T™M = Ker , 

so that the Hessian is non-degenerate on the normal space (TmAT)^. 

This is exactly the non-degeneracy that we need in order to apply the station- 
nary phase Theorem. 

To perform the claimed stationary phase argument, we first take a (small) 
parameter 

6» > 0. 

We use a cutoff in time x(t/ 6 *) with y as in ( 12 . 711 . and evaluate the contribution 


£ C ^ ( 0 ) dt = Ot, ,s {en + 


S{0^* {v)xo{q,p)xi{x,y)PN 



y-qt\ 

J 


dtdxdyd^ dijdqdp. 


When the point {x,y, ^,r], q,p,t) is far from the stationary set M, the integral 
is 0{e°°). Close to the stationary set M, using the fact that the integral carries 
over a compact support, we may use a partition of unity close to M, and on 
each piece we may use straightened coordinates {a, (3) G R^d-i-i-fe ^ such that 

{x,y,^,ri,q,p,t) = 7 ( 0 ,/?), where 7 is a local diffeomorphism, with 
{x,y,^,rj,q,p,t) G M a = 0. 


Using such coordinates, we recover a finite sum of terms of the form 
g(5d+2)/2 dxdyd^dydqdp exp ^^«'(a;,2/,^,r7,g,p,t)^ 

S{0^*{'>i)Pn {t,q,p, ^ X 2 {x,y,^,V,q,p,t) 

= g(5d+2)/2 df 3 exp o 7(a, / 3 )^ 

(^S{.)^*{.)Pn oi{a,(3) X3{a,P), (7.25) 
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where U, Q." are bounded, open subsets, and X 2 i Xa s-re cutoff functions. 
Thanks to the non-degeneracy of the Hessian in the normal direction to 
M, for any (3, we have 


det 


o 7 

Da^ 


(0,/3)^0. 


Hence, by the standard stationary phase Theorem, for any integer J, the above 
integral has the asymptotic expansion to order J 


e^k 2)/2 / exp (-4'o 7(0,/3) 


X 07 x 3 ) (0,/3) (7. 


26) 


^^(fc-5d-2)/2QLj+l 

\ k<2J+d+3 




where the Q 2 j’s are differential operators of order 2j. Now, we anyhow have 


VjeN =0(1). 


On the more, Pn is a polynomial of degree < 4iV in its last argument. This 
implies that the ..) in H7.26|l has at most the size 




Hence, taking J large enough (J > 2N will do), we eventually obtain in 17.2611 . 
using the assumption (H) on the codimension k {k > 5d + 2), 


g(5d+2)/2 I (l-x(l)) (X5(l?e)5„0e(-t)4)dt 

= Oe,Ti .5 (^e^k-5d-2)/2^ ^ 27) 


Fifth Step: Elimination of times such that Tq e < t < 9 - proof of part (ii) of 
proposition ^ 

The previous step leaves us with the task of estimating 

7 ^ (0 U,{-t)(j)e^ dt. 

The idea is to now come back to the semiclassical scale, and write 


37 








This term is expected to be small, provided Tq is large enough. Indeed, the 
propagator exp (—(e^A + n^(x))) acting on </>£ is expected to be close to the 
free propagator exp (—it (e^A + n^(0))) on the time-scale we consider. Hence 
the propagator should have size for large values of time, and the above 

time integral should be 0{Tq ^ 0 as Tq ^ oo. 

We give below a quantitative proof of this rough statement, based on the 
exact computation of the propagator exp (—it (e^A -f n'^{x))) obtained in Theo- 
remlHl The proof given below could easily be replaced by a slightly simpler one, 
upon writing the propagator as a Fourier Integral Operator with real phase. We 
do not detail this aspect, since we anyhow had to use in the previous steps the 
more precise expansion of the propagator given by Theorem]^ this theorem has 
indeed the great advantage to give a representation of the propagator that is 
valid for all times. 

From the second step above (see HT.lbll l. we know 

/ ^ (j) {H^)Se,exp (-it (e‘^A + n‘^{x))) dt 

= OTi,5(e°°)+ f X£"¥ f exp{m{et)/e) 

JTo V"/ dBe-i 

y - Qet 


^ dxdyd^dr]dqdp, (7.29) 


SiO?*iv)xoiq,p)xiix,y)PN t,q,p. 


where we drop the dependence of the phase in {x,y, ^,r], q,p). To estimate 
this term, we now concentrate our attention on the space integral 


^ X 5d 

fe{t) := s = 


/R 6 <J 


exp I i 


'b(et) 


Xoiq,p)xiix,y)PN t,q,p. 


siOriv) 

y - qet 


V~e 


dxdydf^dydqdp. 


We claim we have the following dispersion estimate, uniformly in e, 
|/£(i)| < Ce for some Cg > 0, provided Tg < t < 29/e. 

Assuming (TTiTTll is proved, equation shows that 

p29 


'TqS 


Xeit)x 


Xs {He) Se, Ue{-t)(j)e ) dt 


< CgTo-^- 


(7.30) 


(7.31) 


To —>oo 


(7.32) 


in any dimension d > 3, which is enough for our purposes. It is thus sufficient 
to prove rriTTii . 

We have in mind that the integral defining fe{t) should concentrate 

on the set X = y = q = 0, qt = 0, pt = V, p = f,. Also, the present case should be 
close to the “free” case where the refraction index n'^{x) has frozen coefficients 
at the origin n'^{x) « n^(0). For that reason, we perform in (17.3011 the changes 
of variables 

{x - q)/y/e ^ X, {y- qet) I Vs ^y, q^ Vsq, 

i^p + VsV 9 s(£t, Veq,p) + Vsy- 
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We also put apart the important phase factors in the obtained formula. This 
gives 


feit) = I dqdpdr] exp (^it ^'(p, et, ^/eq, ^/e-q)^ G{q,p, q, et, ^/eq, y^q), 

(7.33) 


/R'ld 


up to introducing the phase 

^{p,et,y^q,y^q) ■= ^ J +„2 (^x{s,y^q,p))^ ds 

Vsp • g - S(£t, Veg, p) ■ X{et, ^q, p) 
et 

!- ^/eq- X{et,eeq,p) 

+ V^V -T-, 

et 

together with the amplitude and compactly supported in p, y/eq) 


G{q,p, q, et, ^/eq, ^/eq) := / dxdyd^ exp ■ (q + x) - iq ■ (y + q)) 

/ a;2 r{et,e£q,p)yy\ 

+ *- 2 - ) 

S{p + (E{et, eeq,p) + e£v)xo {V£q,p) 

XiiV^iq + x),X{£t, V£q,p) + y/ey) Pn {t, V£q,P, v) ■ (7-34) 

Now, the idea is to use the stationary phase formula in the p variable in 1(723, 
where t plays the role of the large parameter. We wish indeed to recognize in 
(177^ a formula of the form 


dp exp 



X smooth(p), 


to recover the claimed decaying factor t jn (723. In other words, we wish to 
get the same dispersive properties as for the free Schrodinger equation. This is 
very much reminiscent of the dispersive effects proved for small times in [Dsf] for 
wave equations with variable coefficients, and relies on the fact that « —p^/2 
as et < 9 is small enough. ^ 

In order to do so, we need to get further informations both on the phase 'I' 
and the amplitude G. 

Firstly, the smooth amplitude G is defined in rr^ . It clearly is compactly 
supported in p and y/eq. Also, the gaussian exp(—a::^/2 + iT(et, y/eq,p)y ■ y/2) 
belongs to the Schwartz space S in the variables x and y (recall indeed 

that Imr(e<) > 0, and et belongs to a compact set), uniformly in the compactly 
supported parameters et, y/eq, and p. From this it follows that the amplitude 
G{q,p,q,et, y/eq, y/eq) belongs to the Schwartz space S (R^'^) in the first and 
third variables q and q, it is <7)1° (R"^) in the second variable p, and these infor¬ 
mations are uniform with respect to the compactly supported parameters et, 
y/eq, together with the (non-compact) parameter y/eq. 

Secondly, the smooth phase dt depends upon the small parameter et S [0, 29], 
together with the two position/velocity variables y/eq and p. All of them belong 
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to a compact set. It also depends upon the variable which is not in a 
compact set. On the more, we have the easy first order expansion in the (small) 
parameter et < 29, 

^{p,et,^/eq,^/eT]) = 

- y + - Veq ' (Veg) - '/s'n ■ {P + 0{9)) + O (6»^) . 

Here the remainder terms 0(9), 0(9^), only depend upon the compactly sup¬ 
ported parameters et < 29 and p, y/eq (they do not depend upon y/erf), and 
they are uniform with respect to these variables. Hence, the stationary points 
of the phase (in the p variable) are given by 

-p - V^77(l + 0(0)) + O (02) = 0. (7.35) 

Finally, there remains to observe that the Hessian of the phase in p is 

yy=-Id + O(0). (7.36) 

Upon taking 0 small enough, all these informations allow us to make use of 
the standard stationary phase estimate in p. More precisely, we write, 

G(q,p, r), et, y/eq, y/evi). (7.37) 

For each given values of q and rj, we analyze the integral over p in USB. If 
y/ep is outside some compact set around the support of G in p, integrations by 
parts in p together with the information 17.3511 . allow to prove that the integral 
over p in rrirzii is bounded, for any integer N, by GN,et ^ for some GN,e > 0 
independent of q and p. Hence the corresponding contribution to fg is bounded 
by GM,et~^ as well. Now, for y/ep in some compact set around the support of 
G in p, we may use the information this, together with the stationary 

phase Theorem with the parameters et, y/eq, y/ep in a compact set, establishes 
that the integral over p in is bounded by Get for some Cg > 0, and Cg 

turns out to be independent of q and p. Hence the corresponding contribuition 
to fg in 17.3711 is bounded by Gt~^G as well. 

All this gives the claimed estimate 

\fe(t)\ < Ggt-‘^// 

The proof of proposition [ 7 | is complete. 

8 Conclusion: Proof of the main Theorem 

We want to prove the convergence 

(w/(t>) 
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when the source S and the test function (j) are Schwartz class. Therefore, one 
needs to prove 


+ 00 

e““=* {Ue{t)Se,(t>e) dt , (j)) as £ ^ 0. 

Proposition n] asserts 

■ i‘2To e / t \ 

e Jo ^ V^J = {w°'"\(t>) 

+ Oto (£°) + o , 

where the notation 0(e°) denotes a term going to zero with e, and Otq (£°) 
emphasizes the fact that the convergence depends a priori on the value of Tq . 
On the other hand Proposition [3 asserts 

^ (1 - X) “ ^- 5 ) 

^o(X)+o(.»). 

Now, for very large times and almost zero energies. Proposition |21 shows, for <5 
small enough, and any k, 

1 /■+°° / \ 

- j e-“'‘ {U,it)xs {He) Se,c^e) dt = 0,A^). 

As for large times and almost zero energies. Proposition 0] shows that, for 5 
small enough, k small enough, and Ti large enough, 

I £ e-“'* (^Ue{t)xs (He) Se, 4 ) dt = 0,,s{e) 

Finally, for moderate times and almost zero energies, one has the following two 
informations. First, for 9 small enough, and uniformly in e, we have 

(ifc) * G) ■" 

= o» (iyi75^) ^ 

Second, for any fixed value of 0 > 0, and Ti, 

^ (1 - x) (££j (ye{t)X5 {He) Se, (j^e) dt 

= Oe.Ti ,s (£°) ■ 

All these informations show our main Theorem, upon conveniently choosing the 
cutoff parameters 9, Tq , Ti (in time), S (in energy), and the exponent k (in 
time). This ends our proof. 
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9 Examples and counterexamples 


9.1 The harmonic oscillator 


Given an appropriate potential V (x), and defining the semi-classical Schrodinger 
operator 

He = + yix), 

our main Theorem proves 


1 

e 



exp ( -i-He ) Se,(l)e 


dt - 



^ exp {—it [—A2^/2 -|- F(0)]) S, (j^ dt. 


(9.1) 


Though we used in many places that our analysis requires a potential of the 
form 


V{x) =-n'^{x) =-nl, + 0{{x) ^), 


it seems interesting to investigate the validity of inu when the potential is 
harmonic 


V{x) = V{0) + J2^^l (9.2) 

1=1 

for some frequencies ujj G M, and a given value V (0) < 0. Such a potential does 
not enter our analysis since it is confining. However, it is easily proved that for 
pairwise rationally independent values of the frequencies Uj, the transversality 
assumption (H) page 123 is true for this potential, whereas in the extreme case 
where all ojj’s are equal, this assumption fails. On the other hand, one may use 
the Mehler formula |Ho| (see for a use of these formulae in the nonlinear 
context) to compute the propagator 


exp I — 


-£^A^/2 + iofxj{2 


d 

1=1 


(9.3) 


n 

1=1 


_ 2iTre sui{ujjt) 


1/2 


exp 


2s sin{LL}jt) 


[[x] + yj) cos{ujjt) - 2xjyj] 


(Here we identified the propagator and its integral kernel). 

Surprinsingly enough, using the Mehler formula to compute the limit on the 
left-hand-side of llOl . we may prove that for rationally independent ujj’s, the 
convergence result mi is locally true in this case, for dimensions d > 4, i.e (lOl 
is true with the upper bounded -l-oo replaced by T, for any value of T > 0. 

We do not give the easy computations leading to this result. The idea is the 
following: at each time kir/ajj (fc G Z), the trajectory of the harmonic oscillator 
shows periodicity in the direction j. However, due to rational independence, at 
times kTrjujj, the trajectory does not show periodicity in any of the d — 1 other 
directions. Hence one gets enough local dispersion from these directions to show 
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that the corresponding contribution to the time integral on the left-hand-side 
of 19.111 is roughly 

/ Al + k-R/ujj)/e \ 

ol / ^ 0 , 

\J{-l+k^/u^i)/e J ^ > 

as long as d — 1 > 2, i.e. d > 4. 

Needless to say, in the extreme case where all w^’s are equal, the result in 
(lOl is false, even locally: in this case, periodicity creates a disastrous accumu¬ 
lation of energy at the origin {all rays periodically hit the origin at times kir/oj, 
k G Z). 

To our mind, this simple example indicates that our main Theorem probably 
holds true for less stringent assumptions on the refraction index. For instance, 
a uniform (in time) version of our transversality assumption is probably enough 
to get the result (without assuming neither decay at infinity of the refraction 
index, nor assuming the non-trapping condition). 

9.2 Examples of flows satisfying the transversality condi¬ 
tion 

We already observed that the harmonic oscillator with rationally independent 
frequencies does satisfy the transversality assumption (H). One actually has 
the value fc = 6 d -I-1 (see (|723)) of the codimension in that case. 

It is also easily verified that the flow of a particle in a constant electric field, 
i.e. the case of a potential 

V{x) = xi, 

does satisfy (H) as well, with A: = 6 d -|- 1. 

Coupling the two flows, it is also verified that the potential 

d 

V{x) = xi +'^uj]x'j/2, 

does satisfy (H) as well, with fc = 6d -I- 1. 

Clearly, these examples are satisfactory, in that we may assume that the 
potential has the above mentioned values close to the origin, and we may trun¬ 
cate outside some neighbourhood of the origin so as to build up a potential that 
satisfies the global assumptions we met in our main Theorem. 
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